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We discuss a relativistic model for heat conduction, building on a convective vari- 
ational approach to multi-fluid systems where the entropy is treated as a distinct 
dynamical entity. We demonstrate how this approach leads to a relativistic ver- 
sion of the Cattaneo equation, encoding the finite thermal relaxation time that 
is required to satisfy causality. We also show that the model naturally includes 
the non-equilibrium Gibbs relation that is a key ingredient in most approaches 
to extended thermodynamics. Focussing on the pure heat conduction problem, we 
compare the variational results to the second-order model developed by Israel and 
Stewart. The comparison shows that, despite the very different philosophies be- 
hind the two approaches, the two models are equivalent at first order deviations 
from thermal equilibrium. Finally, we complete the picture by working out the non- 
relativistic limit of our results, making contact with recent work in that regime. 
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Dissipation 

1. Introduction 

Dissipative fiuid phenomena represent a number of challenges for relativistic physics. 
The main reason for this is the simple fact that the classic Navier-Stokes equations, 
which are not hyperbolic, allow instantaneous signal propagation. This is obviously 
not allowed in a relativistic description, where a model must respect causality in 
order to be considered acceptable. Given the fundamental issues involved, it is not 
surprising that the problem of relativistic heat conduction continues to attract in- 
terest. The issues considered range from questions like "Does a moving body feel 
cold?" (Landsberg, 1967) to issues of stability of models for non-equilibrium ther- 
modynamics (Hiscock & Lindblom, 1983, 1985, 1988; Olson & Hiscock, 1990) and 
whether different descriptions can be distinguished by experiment (Geroch, 1995; 
Lindblom, 1996). 

In the mainstream general relativity community the debate has, to some extent, 
been settled since the late 1970s. The key contribution was the work of Israel and 
Stewart, who developed a model analogous to Grad's 14-moment theory, firmly 
grounded in relativistic kinetic theory (Stewart, 1977; Israel & Stewart, 1979a, b). 
This so-called "second order" theory, which extends the pioneering "first order" 
work of Eckart (1940), has since found applications in a number of contexts. In 
particular, in the last few years there has been a resurgence of interest in the model 
arising from the need to describe highly relativistic plasmas generated in colliders 
like RHIC at Brookhaven and the LHC at CERN (Elze, Rafelski & Turka, 2001; 
Muronga, 2004). However, despite the obvious successes of the Israel- Stewart model 
and its various attractive features there are still dissenting views in the literature, 
see for example Garcia-Colin & Sandoval- Villalbazo (2006); Garcia-Perciante et 
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al (2009). Particular objections concern the complexity of the model and the large 
number of dissipation coefficients that are needed to complete it. This is, however, a 
feature that is shared by all models within the extended thermodynamics paradigm 
(Jou et al, 1993) and it is difficult to see how a simpler model can be developed 
without sacrifice of causality or stability. 

Our discussion is motivated by recent efforts to model the dynamics of super- 
fluid neutron stars. This a problem that requires a general relativistic description 
of "multi-fluid" dynamics, and where thermal and dissipative effects are expected 
to impact on observations. Work in this area is motivated by the effort to detect 
gravitational radiation, and the need to understand various oscillation instabilities 
in rotating neutron stars (Andersson, 2003). These instabilities tend to be counter- 
acted by viscosity and it is obviously important to have a quantitatively accurate 
description of the involved mechanisms. Most studies of this problem have been 
carried out in Newtonian gravity, basically because of a feeling that the relativistic 
problem is too "complicated" . There is, however, a growing body of work building 
towards realistic, fully relativistic, models. The present analysis should be viewed 
in that context. 

Wc are not providing a truly original view of the heat conduction problem. 
Rather, we approach the issue within the multi-fluid paradigm that has been suc- 
cessfully applied to superfluid systems. A key ingredient in this analysis is the 
treatment of the entropy as an additional (masslcss) "fluid" component. This idea 
is obviously phenomenological, but we will show that it provides a model with a 
number of attractive features. Most importantly, the formalism is intuitive and 
readily generalised to more complex settings. 

We build on the convective variational formulation developed by Carter (Carter, 
1976, 1989, 1991), see Andersson & Comer (2007) for an introductory review. In 
addition to the intrinsic elegance of an action principle, an appealing feature of 
the variational approach is that once an equation of state is provided the theory 
provides the relation between the various currents and their conjugate momenta (a 
point that is usually not considered in the context of the heat problem). Another 
key advantage of the variatiational derivation is that incorporating additional fluid 
components is straightforward. 

In Carter's macroscopic model for the heat problem one considers two fluxes, 
one corresponding to the matter flux and one which is associated with the entropy. 
These two currents are the fundamental fields appearing in the matter sector of 
the Einstein- Hilbert action. The Lagrangian of the theory is a relativistic invariant 
and hence it should depend only on covariant combinations of the two fluxes, which 
includes the relative flow between them. This encodes the so-called entrainment 
effect, which tilts the momenta with respect to the currents when two or more 
fluids are coupled, and which turns out to be a crucial feature of the multi-fluid 
approach to relativistic heat conduction. In Carter's original work (Carter, 1988, 
1989), the aim was to keep the model as simple as possible by imposing restric- 
tions on the way that the currents interact. On the grounds of simplicity Carter 
ignored the entrainment. This omission resulted in a truncated model with severe 
stability issues (Olson & Hiscock, 1990). Given this problem, it was suggested that 
Carter's approach does not provide a viable alternative to the Israel-Stewart model. 
This argument is, however, flawed. A detailed comparison of the two formalisms 
(Priou, 1991) shows that, at second order in the deviation from equilibrium, the full 
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(entrained) variational approach is essentially equivalent to the Israel and Stewart 
model. 

An important aspect of the present work is that it extends the recent Newto- 
nian model of thermal dynamics discussed by Andersson & Comer (2010) to general 
relativity. In both cases entrainment has a fundamental impact on the dynamics of 
entropy. In fact, it is an essential ingredient that preserves causality and stability. 
Both discussions lead to a generalization of the Cattaneo equation (Cattaneo, 1948) 
and hence a finite speed of propagation of heat. An equivalent, although not identi- 
cal equation arises in the Israel and Stewart theory. The difference between the two 
approaches is in the underlying thermodynamics. Basically, the Israel and Stewart 
formalism is based on the standard equilibrium Gibbs relation, and therefore the 
thermodynamic quantities take their equilibrium values. In contrast, the covariant 
dynamics of the variational approach leads naturally to an extended Gibbs relation 
(analogous to that in many models of extended irreversible thermodynamics), giv- 
ing the thermodynamic quantities a different meaning. Of course, the theories are 
completely equivalent in the case of thermal equilibrium. 

The paper is structured as follows. Section 2 provides the two-fluid variational 
derivation. The general philosophy that we adopt is that of Carter (1991). We obtain 
the equations of motion by imposing conservative constraints on the variations of 
the Lagrangian density, but we use the explicit freedom in the equations of motion 
to inchide momentum exchange and entropy production while keeping the energy- 
momentum tensor divergence-free. Once the equations of motion are obtained it 
is useful to make a choice of frame to discuss the thermodynamics and the Gibbs 
relation. We discuss this point in detail, arguing why the choice that we make is 
natural. Section 3 contains the derivation of the relativistic generalization of the 
Cattaneo equation that follows once we impose the second law of thermodynamics. 
The result is then compared to the predictions of the Israel and Stewart model. 
The Newtonian limit of the theory is obtained in Section 4, which establishes the 
close connection to the non-relativistic model of Andersson & Comer (2010). The 
paper concludes with a brief discussion of the implications of the results, possible 
future extensions and applications. 

Throughout the paper we use the convention that 4-dimensional spacetime in- 
dices are represented by lowercase letters starting from the beginning of the al- 
phabet, a, h, c, ... while 3-dimensional spatial indices are lowercase letters k, .... 
There should be little risk of confusion. We denote the covariant derivative by a 
semi-colon. 

2. The two-fluid model 

We consider the problem of heat conduction in general relativity at the macroscopic 
level. This means that we assume that the particle number is large enough that 
the fluid approximation applies and that there is a well defined matter current, n". 
Moreover, we adopt the multi-fluid formalism developed by Carter (1989) and treat 
the entropy as an effective fluid with flux s". This current is in general not aligned 
with the particle flux. The misalignment is associated with the heat flux and leads 
to entropy production. 

For a generic two-fluid system, the starting point is the deflnition of a rela- 
tivistic invariant Lagrangian-type master function A. Assuming that the system is 
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isotropic, we take A to be a function of the different scalars that can be formed by 
the two fliixesf. Prom n" and s" we can form three scalars; 

(2.1) 
(2.2) 
(2.3) 

Hence, we can write the master function as the density for the matter sector of the 
Einstein-Hilbert action; 
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An unconstrained variation of A then leads to 
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(2.4) 
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Using (2.1)-(2.3), we can change the passive density variations for dynamical 
variations of the worldlines generated by the fluxes and the metric (Andersson & 
Comer, 2007). That is, we use 
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This means that the variation (2.5) becomes 
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Prom (2.9) we can read off the conjugate momentum associated with each of 
the fluxes; 



Ma 



dA 
dn 

dA 



gabiB's'' + A'^'n') , 



where we have introduced the coefficients 



an^ as^ dj'' 



(2.10) 
(2.11) 

(2.12) 



t It should be noted that we consider the simplest "convective" model. The natural way to 
a<;count for viscosity would be to allow the master function to depend also on the associated 
stresses, see Carter (1991) and Prion (1991). This model is, however, significantly more complex 
and we do not consider it here in order to keep the discussion clear. 
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The conjugate variables (2.10) and (2.11) demonstrate the fundamental role 
of the master function (Carter, 1989). The distinct roles of the fluxes and their 
conjugate momenta are often not considered in the fluids literature. A key advantage 
of the variational approach is that the quantities are immediately determined by 
the form of the master function. Moreover, it is clear that the momenta are not 
generally aligned with the respective currents, due to the fact that the master 
function depends on the relative flux. This is an important effect. Fundamentally, 
there is no physical argument to rule out the dependence on j^. In fact, this coupling 
is associated with the cntrainment effect that is known to be of central important 
in other multi-fluid systems. 

In the case of superfluid neutron star cores the entraiment arises due to the 
strong interaction and couples the neutron and proton fluxes, sec for example Comer 
& Joynt (1988). In the present problem, with matter and entropy, we will show that 
the entrainment is associated with the thermal relax;ation of the system. This is an 
important effect that must be accounted for. 

The energy-momentum tensor is obtained from (2.4) by noting that the dis- 
placements of the conserved currents induce a variation in the spacetime metric 
and therefore the variations of the fluxes, (5n" and Ss"', are constrained, sec Carter 
& Quintana (1972); Andersson & Comer (2007); Prix (2004) for discussion. The 
energy-momentum is thus found to be 

T^'> = Han'> + eas'' + m„^ , (2.13) 

where we define the generalized pressure, ^, as 

* = A - M„n" - ^asV (2.14) 

The equations of motion are obtained by requiring that the divergence of the 
energy- momentum tensor (2.13) vanishes. For an isolated system, we can express 
this requirement as an equation of force balance 

Ta'.b = /a + /: = , (2.15) 

where the individual force densities are 

/: = 2^[„;6]S^ + sV„ . (2.17) 

The square brackets represent anti-symmetrization in the usual way. 

We note that, in order to obtain the energy momentum tensor (2.13) we needed 
to impose the conservation of the fluxes as constraints on the variation. However, 
the equations of motion, (2.16) and (2.17), still allow for non-vanishing production 
terms. If wc, for simplicity, consider a single particle species, the matter current is 
conserved and we have 

<a = . (2.18) 
This removes the second term from the right-hand side of (2.16). In contrast, the 
entropy flux is generally not conserved. Wc will have 

s% - P. > , (2.19) 

in accordance with the second law of thermodynamics. This suggests that, to make 
progress, we need to connect the variational results with the relevant thermody- 
namical concepts. In doing this it makes sense to make a specific choice of frame. 



Article submitted to Royal Society 



6 



C.S. Lopez- Monsalvo and N. Andersson 



(a) The m,atter frame 

The conservation law (2.18) implies that the force is orthogonal to the matter 
flux, n", and therefore has only three degrees of freedom. Furthermore, because of 
the force balance (2.15), wc also have n°/^ = 0. This suggests that it is natural to 
focus on observers moving with the matter frame. We associate the matter current 
with a four-velocity such that 

n° = mi", (2.20) 

where UaU°- = — 1 and n is the number density measured in this frame. Historically, 
this is known as the Eckart frame. As will soon become clear, this is the natural 
frame in the case of a single particle species, essentially because it simplifies the 
analysis. More complex settings, e.g. when dealing with additional particle species 
and reactions, will make the choice of frame less obvious. It may well be that the 
best strategy in such cases is to follow Landau & Lifshitz (1959) and work in the 
centre of mass frame. Having said that, it is worth noting that even in the more 
general problem is there a unique frame associated with the entropy/heat flow. 
Even though we will not discuss this problem further here, there are interesting 
issues that warrant more detailed thinking. 

Having chosen to work in the matter frame, we can decompose the entropy 
current and its conjugate momentum into parallel and orthogonal components. 
The entropy flux is then expressed as 

s« = s*(n" + «;°) (2.21) 

where is the relative velocity between the two fluid frames, and u'^Wa = 0. 
Letting = su^ where is the four- velocity associated with the entropy flux, we 
see that s* = s'y where 

7= + = (l-W;2)-l/2 ^ (2.22) 

is the redshift associated with the relative motion of the two frames. In the following, 
we will use an asterisk to denote matter frame quantities. 
Similarly, we can write the thermal momentum as 

Oa = e*Ua + e^Wa = (B^S* + A'^^Tl) Ua + B's*Wa , (2.23) 

where we have made use of (2.11). From these expressions we readily obtain a 
measure of the temperature measured in the matter frame; 

-u"ea = e* = B^s* + A'^^n . (2.24) 

In essense, this measure respresents the effective mass associated with the entropy 
component. We have also deflned 

= B^s* . (2.25) 

It is worth noting that, if we ignore the coupling between the fluxes in the mas- 
ter function by taking A'^^ = 0, then we have 6* = ^'f. This particular case was 
considered by Carter (1988), and from the analysis of Olson & Hiscock (1990) we 
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know that it leads to a model that exhibits instabilities. It is useful to keep this 
in mind during the following developments. As we will see, the main problem with 
Carter's "regular" model is that it leaves no freedom to adjust the thermal relax- 
ation timescale. 

In order to express the energy-momentum tensor in terms of the matter frame 

quantities, wc define the variables a"" = s*w°' and pa = B^s*Wa. Using the above 
expressions the stress-energy tensor (2.13) can be written in a more familiar form; 

Tab = - [A - UaUb + 2u(„g6) + P„(, , (2.26) 

where, making use of the projection orthogonal to the matter flux 

hab = 9ab + UaUb , (2.27) 

the heat flux (energy flow relative to the matter) is given by 

Qa = -habUcT''^ = S*e*Wa • (2.28) 

We also have 

Pab = hab'if + PaCTb . (2.29) 

It is worth noting at this point that the variational analysis leads to the presence of 
"shear" terms in the energy- momentum tensor. Such terms arc usually associated 
with viscous stresses, and it is interesting to note that they arise even though we 
consider the pure heat conduction problem. Moreover, this exercise shows that the 
energy density measured in the matter frame can bo obtained by a Legendre-type 
transform on the master function. That is, we have 

p* = UaUbT"^ = -A + p^a" . (2.30) 

In fact, this relation informs the choice of cr" and Pa as key variables (Carter, 1976). 

(6) The temperature problem 

Thermodynamic properties such as pressure and temperature are uniquely de- 
fined only in equilibrium. Intuitively this makes sense since, in order to carry out 
a measurement (of say the temperature), the measuring device must have time to 
reach "equilibrium" with the system. The measurement is obviously only meaning- 
ful as long as the timescale required to obtain a result is shorter than the evolution 
time for the system. Of course, this does not prevent a generalisation of the various 
thermodynamic concepts. The procedure may not be unique, but one should at 
least require the generalised concepts to be internally consistent within the chosen 
extended thermodynamics model. As a useful demonstration of this notion, and the 
fact that our model satisfies this criterion, we will consider the particular case of 
the temperature. 

When the system is out of equilibrium, we can define a number of different 
"temperatures" . It makes sense to refer to the quantity obtained from (2.24) as the 
dynamical temperature since it corresponds to the effective mass of the entropy 
component. Similarly, the effective mass for the particles, /i, is given by 

H = -u"fia = B'^n + A'^^s* . (2.31) 
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Let us now show that (2.24) agrees with the thermodynamical temperature that an 
observer moving with the matter would measure. Using the standard definition of 
temperature, we consider the variation of the energy with respect to the entropy in 
the observer's frame (while keeping the other thermodynamic variables fixed). To 
do this, we note that the energy density measured in the matter frame (2.30) is a 
function of three independent state variables, p* = p* (n, s* ,p). 
Determining the energy density directly from (2.13), we get 

p* = ^n + e*s* . (2.32) 

Using the definitions (2.24), (2.25) and (2.31) we can evaluate the generalized pres- 
sure (2.14) in the matter frame as a Legendre-type transform of the master function; 

'^{p,0*,p) = A + nn + 0*s* -pa , (2.33) 

and so the variation of (2.32), 

dp* = /idn + 0*ds* + adp , (2.34) 

shows that the dynamical temperature agrees with the thermodynamical tempera- 
ture provided we evaluate it in the appropriate frame. We also see that, when the 
system is out of equilibrium, the energy variation (2.34) depends on the heat flux 
(encoded in (7 i uid p„). This extended Gibbs relation is similar to that which is 
postulated in many approaches to extended thermodynamics (Jou et al, 1993). The 
main difference here is that (2.34) arises naturally from the variational analysis. 

This result is far from trivial. The requirement that the two temperature mea- 
sures agree determines the additional state parameter, p, to be held constant in the 
variation of p*. Any other choice of the third parameter, e.g. j^, will lead to the 
determined temperatures being different and, hence, the model less consistent. Ad- 
ditional evidence that we have identified natural parameters of the non-equilibrium 
problem is provided by the Legendre transformations (2.30) and (2.33). This im- 
portant point was originally made by Carter (1976) in a work that is not widely 
known. 



(c) Thermal equilibrium 

Since the main part of our discussion concerns the dynamics of systems out of 
equilibrium, and the comparison of different possible models, it makes sense to make 
a few comments on the state of equilibrium. As usual, we take thermal equilibrium 
to mean that there is no heat flux. Hence, we have q"^ — and the entropy is carried 
along with the matter. It then makes sense to introduce the equilibrium quantities 
p, s and 6, evaluated in the matter frame, in terms of which we recover the standard 
Gibbs relation; 

dp = pdn + eds . (2.35) 
In equilibrium, both n" and are conserved which means that 

<a = — = — , (2.36) 

' n s 

where 

n = u^n-a , (2.37) 
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and similar for s. 

It is worth noting that if we introduce S = s/n then we have 

5 = 0. (2.38) 

That is, the specific entropy remains constant in the matter frame. 
By adding the momentum equations (2.16) (2.17) we find that 

{ufi + se)u + h\ {ni2-b + sO-fi) = . (2.39) 

Making use of the fundamental relation 

P + p = niJ, + se , (2.40) 

where P is the equilibrium pressure (to be distinguished from the generalized pres- 
sure that is relevant also out of equilibrium), and 

dP = ndfj, + sde , (2.41) 

we have 

hi [Pb + {P + p)ub] =0 {P + p)u'' = -hlPb . (2.42) 

As expected, we have the usual relation between the acceleration and the pressure 

gradient. Finally, the entropy momentum equation can be cast in the form (the 
easiest way to see this is to set Fg = in the analysis that follows in the next 
section) 

h"'' {VbO + 0Ub) = -^611" = -hle-b ■ (2.43) 
Comparing these last two expressions we see that we must have 

In other words, temperature changes lead to pressure variations and vice versa. 

3. The relativistic Cattaneo equation 

It is well-known that the classic thermodynamics description leads to non-causal 
heat conduction. This is obvious if we consider the general solution to the heat 
equation corresponding to given initial data in an unbounded region of space. Due 
to the parabolic nature of the heat equation, any initial profile evolves to predict a 
non-zero temperature throughout space, even at arbitrarily early times. To remedy 
this problem, and impose a maximal propagation speed for heat, Cattaneo (1948) 
proposed a modification to Fourier's law such that the heat flux vector, q, is related 
to the temperature gradient according to 

rq-hq=-fcVT. (3.1) 

This model incorporates the (finite) thermal relaxation time t expected to arise 
from the fact that, on the micro-physical scale, heat propagates due to particle 
collisions. Cattaneo's equation leads to the temperature satisfying a telegraph-type 
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equation whose hyperbolic nature leads to a finite propagation speed of heat pulses 
(Jou et al, 1993). 

Since any rclativistic model must encode causal heat propagation, one would 
expect the heat flux to be described by an equation similar to (3.1). The antici- 
pated form for the relativistic analogue of Cattaneo's equation is, indeed, generally 
agreed upon but its derivation and the physical interpretation of the involved vari- 
ables differ among proposed models. Our main aim is to obtain the relevant heat 
conduction equation within the variational approach, and compare the result to the 
second order model of Israel & Stewart (1979b). In doing this, we will highlight the 
differences between the two models. 



(a) The variational approach 

Wc want to formulate; a rclativistic analogue of Cattaneo's equation. The basic 
strategy will be to use the orthogonality of the entropy force density with the 
matter flux, solve for the entropy production rate Fg and flnally impose the second 
law of thermodynamics. 

Let us first note that, by making use of (2.28) we can express the entropy flux 
in terms of the matter flow and the heat flux; 

= s*"" + ^9" ■ (3-2) 
Moreover, the conjugate momentum becomes 

ea = 6*Ua+l3qa , (3.3) 

where /3 is given by 

n 9^-9* \ / 1 .A^^nX 

\s* ~^ s*e* ) ~ \s* s*e* ) ' ^ ' ' 

It is worth noting that if wc set A^^ = 0, as in Carter's "regular" model (Carter, 
1988), then /3 takes the specific value 1/s*. 

In terms of the variables we have introduced we have 

rs», = 20[a;6]«"a'^ , (3.5) 

which leads to the entropy creation rate 

r, = S% = -^q" {e*a + 6*Ua- Pqc;aU^ + Pqa + AcW'^a) • (3.6) 

Dots represent time derivatives in the matter frame, as before. That is, is the 
four-acceleration and 

- u'^q'',, , (3.7) 

represents the time variation of the heat flux. 

The expression given in (3.6) has to satisfy the second law of thermodynamics. 
The simplest way to achieve this is to demand that the entropy production is a 
quadratic in the sources. This suggest that the heat flux takes the form 

= -nh"'' {e*.t, + e*Ub + 2/3g[b.,]U^ + $qb) (3.8) 
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where k > is the thermal conductivity of the fluid. This leads to a relativistic 
analogue to Cattaneo's equation, which can be written 

2r9[»:=l«c + 9" = + ' (3-9) 

or 

T (g" + gew^'«) + g» = -K/i"*" + r«6) , (3.10) 

where we have used the fact that Uc?"'" = —qcU"''". Moreover, we have introduced 
the effective thermal conductivity k; 

K=—!^, (3.11) 

and the thermal relaxation time 

r = ^ . (3.12) 

It is worth noting that, if /3 varies on a timescale (say) which is long compared 
to the thermal relaxation, then 

K/5 « ^ » ^ < 1 . (3.13) 

One would expect this to be the case in most situations of practical interest. Thus, 
we would simply have 

KKi K , and T !v kI3 . (3-14) 

(&) A different perspective 

Before we proceed to compare (3.9) to the corresponding equation obtained 
within the Israel & Stewart (1979b) approach , it is worth considering the varia- 
tional result from a slightly different point of view. As in the Newtonian problem 
(Andersson & Comer, 2010), the entropy momentum equation provides an evolu- 
tion equation for the heat flux. This equation can be written in a number of ways, 
some of which are helpful in interpreting the results. 

As a first step, we recall that is orthogonal to m". By contracting (2.16) with 
and using the result in (2.17) we arrive at an elegant expression for the force 

-e*r, = 2u^St + (3.15) 

Written in this form, the force is clearly orthogonal to w". Moreover, this expression 
emphasises the relevance of the entropy momentum 9a- However, if we want to gain 
insight into the key factors that contribute to the force then we need to expand this 
expression. To do this, we contract equation (2.17) with (3.2) to get 

(.''W = ^V^ = ^9V^ (3.16) 

This shows that the entropy production only depends on the piece of the force 
which is parallel to the heat flux. In general, we can decompose the entropy force 
into two terms; 

fa = fUa + fi, (3.17) 
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where both pieces are orthogonal to m", and is also orthogonal to g". From 
(3.16) it is obvious that cannot contribute to the entropy production. Hence, 
this term is not constrained by the second law. This is an important point because 
there is no obvious way to distinguish the viability of models with different forms 
of . Given this, it is interesting to consider the specific force terms that arise in 
the variational formalism. 

With the definitions above it is straightforward to show that 



^ /Ilk i /I* 



Meanwhile, from (3.16) we have 



where 



r =J-^f\\ 



S^Oa = S*0* + _ 



(3.18) 



(3.19) 



(3.20) 



As expected, this takes us back to (3.6). Moving on to the non-entropy producing 
part of the force, the Ansatz (3.8) together with a projection orthogonal to leads 
to 



= 

J a 



where 



LC 6 



- Wqc).t 
2 



s*6* 



PqUb-.d 



qcq 

/y2 



, (3.21) 



(3.22) 



projects out the component that is orthogonal to We see that the variational 
approach leads to the presence of terms that, even though they involve the heat flux, 
are not associated with entropy production. As far as we are aware, the dynamical 
role of these terms has not been discussed in detail in the literature even though 
similar terms are (as we will soon see) also included in the Israel-Stewart formalism. 
The variational model leads to these terms taking a specific form. In particular, the 
terms in (3.21) are all quadratic in q°-, the deviation from equilibrium. At this order, 
the most general case would allow a force of form 



(3.23) 



with Aa, B and C unspecified coefficients. There may also, in principle, be first 
order terms. Clearly, (3.21) represents a particular case where all the coefficients 
follow from (3. Hence, the variational model is a particular example of the general 
class of permissible theories. The fact that all these models satisfy the second law 
of thermodynamics means that we cannot express a preference at this point. A 
very interesting question concerns whether there are situations where has a 
distinguishable effect on the dynamics of the system. If one could show that this is 
the case, then we may be able to narrow down the possibilities. 
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(c) Remarks on the role of the four acceleration 

It is worth commenting on the presence of the term associated with the four- 
acceleration on the right-hand side of (3.10). As we will see later, this term has no 
counterpart in the Newtonian problem. However, its presence in the rclativistic heat 
equation has been known since the pioneering work of Eckart (1940). Formally, this 
term originates from the local energy balance, (3.31). Physically, it results from 
the fact that the infinitesimal 3-spaces orthogonal to the matter world lines are 
not parallel, but relatively tipped over because of the curvature of the world line. 
This leads to the interpretation of the four-acceleration contribution in terms of 
the effective inertia of heat (Ehlcrs, 1973). Interestingly, Garcia-Colin & Sandoval- 
Villalbazo (2006); Garcia-Perciante et al (2009) have recently suggested that this 
term may be the origin of instabilites. We will not discuss this suggestion in de- 
tail, but note that there are varying points of view in this area of research. Most 
researchers seem to accept that the four acceleration term is both inevitable and 
physically meaningful. 

The derivation of (3.9) was based on the thermal momentum equation (2.17). 
We could, in principle, also make use of the other momentum equation, (2.16). 
Contracting this equation with s" we have 



where 



. 2a . 

IJ-Ua + li-,a + qaOl H U gu-w 

n ^ ^ 



(3.24) 



a = = -T— . (3.25) 

This gives us an expression for q°"Ua which could be used in (3.18). We would then 

arrive at a different form for Fg and as a result, the Cattanco equation will also be 
different. The four acceleration term will now be replaced by the chemical potential 
gradient. After some algebra, we arrive at 



X q" + jH;b + + joj + 2 - joj u\,,b] I • (3.26) 

We could impose the second law on this result, and hence derive an alternative 

form for the Cattaneo equation. There is, however, no obvious reason why this 
form would be more useful than (3.9). Hence, we will not pursue this strategy 
further here. 



(d) The Israel-Stewart approach 

The most successful approach to the problem of causal heat conductivity and 
dissipation in relativistic fluid dynamics is due to Israel and Stewart (Stewart, 1977; 
Israel & Stewart, 1979a, b). A detailed comparison between their results and the 
variational model has already been carried out by Prion (1991). The key results of 
this comparion are: i) The two models differ only at second order in the deviation 
from equilibrium, ii) The inclusion of entrainment, A'^^ ^ 0, is essential in the 
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variational analysis, iii) The two models belong to a larger class on non-equilibrium 
thermodynamics models. We will demonstrate these results in the particular case 
of pure heat conductivity. Prion's analysis includes the viscosity contributions from 
Carter (1991), which makes many of the equations rather complex. The message is 
clearer if we focus on the heat conductivity problem, and we feel that it is important 
to make the comparison as transparent as possible. 

In order to effect the comparison, it makes sense to begin by working through 
the derivation of the Israel-Stewart model for a heat conducting system. We focus 
on the phcnomcnological description, and simply note that the model has a firm 
foundation in relativistic kinetic theory. Hence, we take as the starting point the 
stress-energy tensor (Hiscock & Lindblom, 1983); 

= pu^u^ + {P + T)h''^ + 2u("g*) + (3.27) 

The main difference from the variational model is that the thermodynamical quan- 
tities refer to an equilibrium state. As before, p is the energy density, P is the 
pressure and is the heat flow orthogonal to the matter flow. Meanwhile, t"''' 
and T are the stresses caused by viscosity in the fluid. The tensor r"'' satisfles the 
relations 

= = = T[„,] . (3.28) 

Motivated by kinetic theory, Israel and Stewart expands the entropy flux to 
include a complete set of second order terms 

„o 1 „.a a t"" 

s- = sn-+'L--- (p,r^ + jj^q\, + /J^roer"^) _ + + . (3.29) 

Here, T is the absolute temperature associated with the equilibrium state and the 
coeflacients /?o, /32, o-o and ai correspond to different couplings that need to 
be provided (i.e. obtained from the microphysics of the problem). By neglecting 
all viscosity contributions we are left only with the /3i term, and the entropy flux 
(3.29) reduces to 

= sn" + Ig-^ - ^q'u\ (3.30) 

From this definition it is straightforward to calculate the divergence of sf to 
impose the second law. Prom the equations of motion we also have 

u,T'^\, = -p - i(p + P)u-.^ - q-,^ + = . (3.31) 

We now use the fundamental relation of thermodynamics in the equilibrium 
form 

P + p = pn-\- sT , (3.32) 
together with the equilibrium equation of state, p = p{n, s), which leads to 

P;a = l^n.,a + Ts,a • (3.33) 

We also need to use the fact that the flux n" is conserved. 
Combining these results, we arrive that 



,1 L^b 

5/;a — 



T,b + Tiib + TPiqb + [ i^u" ] T^qb 



2T 



;a 



(3.34) 
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As before, we need to ensure that the entropy production is positive. This follows 
as long as 



q" = -kTh 



ab 



T-b + Ub + Mb + 



T 



T 

— L 

2 



(3.35) 



This relation provides the Israel-Stewart version of the relativistic Cattaneo equa- 
tion. 

As in the variational case, one may add other terms to the heat flux as long as 
they do not l(;ad to the generation of entropy. As an example, Hiscock & Lindblom 
(1983) included the term 

-72ftT/i»''(7^«[e;6] , (3.36) 

the presence of which is motivated by kinetic theory. This term is clearly orthogonal 
to both and and, hence, does not affect Fg. Technically, it is also a second 
order term because the fluid shear is associated with the viscosity. 



(e) The Extended Irreversible Thermodynamics view 

We have reached the point where we can compare the final equations for the 
heat flux from the variational approach, Eq. (3.9), to the Israel-Stewart model, 
Eq. (3.35). However, before we do this, let us consider the problem from the "ex- 
tended irreversible thermodynamics" point of view. This approach, which was first 
developed by Carter (1976), provides an immediate illustration of the fundamental 
difference between the two models that we are discussing. 

As in the variational analysis, we take the heat flux to be given by 

s- = s*u'' + . (3.37) 
The corresponding stress-energy tensor is 

T"'' = p*u''v'' + 2?i("q'') + + cr V , (3.38) 

where \E' denotes the (generalised) pressure, and the variation of the energy density 
(in the matter frame) is 

d/>* = /idn -I- e*d.s* + adp . (3.39) 

Obviously, we are now refering to out-of-equilibrium quantities. Hence, we use an 
extended Gibbs relation that accounts for the heat flux. This immediately leads to 

p% = /xn;a + + ^(/9g);a . (3.40) 

Using these relations and the particle conservation, we find that 

UbT% = {n^i -p* - *)e -6*s* -ap + qb-^aV^u" - q% + pK^^a^Ub . (3.41) 

At this point we need the fundamental relation (2.32), which allows us to complete 
the derivation. We arrive at 

s% = -^Q''[(^*;b + 6*<^b-qbP + {Pqb);cu'] , (3.42) 

which is identical to (3.6). Hence, the heat flux will (again) be governed by (3.9). 
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(/) Comparing the results 



Since the two approaches are based on different strategies, any comparison be- 
tween the variational model and the Israel-Stewart results must be done carefully. 
Notably, the Israel- Stewart model is based on an expansion including terms up to 
second order in the deviation from equilibrium. Meanwhile, the variational analysis 
did not involve such an expansion. As a result, the final equation for the heat flux, 
(3.9), contains higher order terms while (3.35) is manifestly linear in (j". 

In order to compare the results, it makes sense to focus on the linear deviation 
from equilibrium. Then we have 



r = r + o(g2) , 

and it follows that (3.10) can be approximated by 



(3.43) 



(3.44) 



Here it is worth noting that in the second term in the first bracket we could use 
the standard decomposition of Ua-.b in terms of the shear, expansion etcetera. This 
would lead to terms that were explicity excluded from the Israel-Stewart model at 
the point where we focussed on the case with Tab = r = 0. Basically, the variational 
analysis leads to the presence of terms that couple the heat flux to the shear and 
expansion of the flow. As these were artificially excluded from the analysis that 
led to (3.35) we cannot count this as a difference between the two models. In fact, 
the full comparison carried out by Prion show that these terms agree in the two 
descriptions. 

Keeping terms up to second order (treating the shear and the divergence of 
as first-order quantities), (3.35) can be written 



r/g" + g« « -K7/i«* {T,b + Tub) 



where 



and 



2 



1 + m" 



kT^ 

2 



(3.45) 



(3.46) 



(3.47) 



Now it is clear that the two equations for the heat flux are formally identical, 
and we can "identify" the parameters in the two models. The upshot of this is 
that, the models will only produce different results at higher order deviations from 
equilibrium. Given that this regime is hardly tested at all, we cannot at this stage 
comment on which of the two descriptions (if either) may be the most appropriate. 
Having said that, it is clear that the variational approach is formally elegant and 
the fact that it applies also far from equilibrium (at least in principle) may be 
relevant. An interesting question concerns whether there are situations where the, 
rather specific, set of higher order terms predicted by the variational analysis affect 
the nonlinear dynamics. 



Article submitted to Royal Society 



Thermal Dynamics in General Relativity 



17 



Before we conclude the comparison, it is worth noting that the difference be- 
tween the two models was apparent already from the beginning. This is clear if we 
consider the stress-energy tensor. In the variational case we have 

Tab = P*UaUb + 2q(aUb) + "^Kb + PaC^a ■ (3.48) 

Comparing this to (3.27) we see that 

T=\{s*fB''w'' (3.49) 

and 

Tab = 3(s*)2r {waWb - Kbw'') (3.50) 

These terms are quadratic in w°' (that is, q°'). Hence, it is obvious that, in order to 
carry out detailed a comparison we ought to include also shear- and bulk viscosity 
in the model, c.f. the analysis of Prion (1991). 



4. The Newtonian limit 

Having developed a model for heat conductivity in general relativity, and discussed 
the results in the context of the well-established Israel-Stewart model, we should 
also consider the implications of the model for non-rclativistic systems. The close 
connection between the variational multi-fluid approach in Newtonian gravity and 
extended irreversible thermodynamics has already been discussed by Andersson & 
Comer (2010). The results demonstrated that a two-fluid model based on a mas- 
sive component representing the particles and a massless component representing 
the entropy, reproduces many key results from the literature (Jou et al, 1993). In 
particular, the non-relativistic Cattaneo equation is obtained immediately from the 
momentum conservation law for the entropy. The non-relativistic analysis is, in 
fact, completely analogous to the discussion in this paper. To demonstrate this, 
and relate the two models, we will now work out the Newtonian, low- velocity, weak 
gravity, limit of our main equations. 

Let us first return to the variation of the master function A. Assuming low 
velocities we have 



f = -n^Sa « sn 1 -h 



w 

2? 



(4.1) 



where = w^^ represents the (squared) magnitude of the relative velocity between 
matter and entropy (w^g = — vl). Note that we need to keep the speed of light, 
c, explicit in this discussion. This then leads to 



s \ 1 + ] dn + n 



(4.2) 



and 
dAR 



dn - 



ds-^A'^'dw^ . 

(4.3) 

To make progress it is essential to appreciate that the Newtonian limit is singu- 
lar, see for example the rigorous analysis of Carter & Chamel (2004). In order to 
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effect a direct "calculation" it is useful to separate the "ballistic" rest-mass contri- 
bution to the master function. That is, we use (recalling that the entropy is taken 
to be massless) 

A = —mnc^ — E{n, s, w^) . 
From the above results it follows that 



OA 
dn 



nS" + s ( 1 + ^ M" 



defining the chemical potential, /j., and 



dA 



w 

2? 



-mc — 



= T 



(4.4) 
(4.5) 



(4.6) 



which defines the temperature T (as the entropy chemical potential). Finally, 



dA 



26 



dE 



a 



(4.7) 



defines the Newtonian entropy entrainment parameter a (Prix, 2004). These three 
relations allow us to express B", and A'^^ in terms of the Newtonian coefficients 
/U, T and a. 

We also need the weak field spacetime metric. To the required order, we have 
the line element 



1 + 



2$ 



dt^ + Qijdx^dx-' 



(4.8) 



where $ is the gravitational potential and gij is the flat space metric. Meanwhile, 
the four velocities have components 



1 



and 



-^ 


h — 




2c2 








2c2 



u„ = v„ 



u» = V, 



This means that the two fluxes become 



nu: 



n = 




_ n 



1 

nc 



2c2 



(4.9) 

(4.10) 

(4.11) 



Using the different expressions in the momentum equations, keeping only terms 
of order unity, we obtain the Newtonian equations of motion. For the particles we 
then flnd 

/ 2a 
nV?.„ « nidt + viVj) mv^ + —w 



-nVi(At + TO$) - 2awi^ViVj , 
This means that the relevant momentum equation becomes 

n{dt + viWj) (mvl + —wf] + nS7i{n + m$) + 2awi^Wiv'^ = 



(4.12) 
(4.13) 

(4.14) 
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For the massless entropy, it is easy to see that we get 

s-f,l.^^s{dt+viVj)(^^wr^ , (4.15) 

syi,, « -sViT - 2awi,Viv'j , (4.16) 
and the final momentum equation becomes 

s{dt + viVj) + sViT + 2awi,Wiv'j = . (4.17) 

In order to facihtate a direct comparison with the discussion by Andersson & Comer 
(2010) we now use the entropy conservation law 

s''.^^dtS + Vj{svi)=T, , (4.18) 

and define 

7rf = 2a<" . (4.19) 

Then (4.17) can be written 

dt^i + [viT^!) + sViT + Tv^^Vivi = ^nt = f! , (4.20) 



or 



(5*+<)7r|- V, ( ^ j -TT^V, i^j+sViT + n^jVivi + n^^V^vi = ^tt? . (4.21) 



We have arrived at the Newtonian two-fluid point model that was taken as 
the starting point by Andersson & Comer (2010). In other words, their Newtonian 
model is the natural non-relativistic counterpart to the model developed in Sec- 
tion 2. That this had to be the case was, more or less, obvious given the discussion 
by Andersson & Comer (2007), but it is still useful to have a direct comparison. 
After all, the derivation shows explicitly that the four acceleration term in (3.9) is 
a purely rclativistic effect. The comparison also clarifies the; physical interpretation 
of the different variables, and the meaning of the various parameters. 



5. Discussion 

We have discussed a relativistic model for heat conduction. The model builds on the 
convcctive variational approach to multi-fluid systems designed by Carter (1989), 
and focusses on the role of the entropy as a dynamic entity. The model assumes 
that the entropy can be treated as a "fluid" distinct from the matter flow. We have 
demonstrated how this approach leads to a rclativistic version of the Cattaneo equa- 
tion, encoding the thermal relaxation time that is needed to satisfy causality. More- 
over, we have shown that the model naturally includes the non-equilibrium Gibbs 
relation that is a key ingredient in most approaches to extended thermodynamics. 
By focussing on the pure heat conduction problem, neglecting other "dissipation 
channels" , we compared the variational results to the celebrated second-order model 
developed by Israel & Stewart (1979b). The comparison showed that, despite the 
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very different philosophies behind the two approaches, the two models are equiv- 
alent at first order deviations from thermal equilibrium. This was not surprising. 
In fact, Priou (1991) has already carried out a similar comparison of the corre- 
sponding second-order models (including viscosity). His results show that the two 
models contain the same key elements, and that they belong to a wider class of 
permissible models. The simpler context of our analysis serves to clarify the main 
points. Finally, we worked out the non-relativistic limit of our results, making con- 
tact with the recent work of Andersson & Comer (2010). This essentially completes 
the picture, and we now have a consistent framework for discussing causal heat 
conductivity in both Newtonian and relativistic dynamics. 

Our discussion obviously only scratched the surface of what is a very difficult 
problem. Wc did not address foundational issues concerning the link between this 
kind of phenomenological model (e.g. the entropy "fluid") and the relevant micro- 
physics/statistical physics. This is a rich and challenging area, where many issues 
remain to be resolved, and one can imagine several interesting directions in which 
the current work may be developed. This work may also be applied in a number 
of exciting contexts, ranging from high-energy collisions probed by, for example, 
RHIC and the LHC, to neutron star dynamics and issues relevant for miiltimessen- 
gcr astronomy, and even cosmology and the evolution of the Universe itself. 

As this work was carried out, three truly inspirational colleagues passed away. Without 
the work done by these pioneers our understanding of relativistic heat problems would 
not be what it is today, and we would like to dedicate this work to the memories of Peter 
Landsberg, Bill Hiscock and Jiirgen Ehlers. 

NA acknowledges support from STFC via grant number ST/H002359/1. CSLM grate- 
fully acknowledges support from CONACyT. 
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